Small Angle Muon and Bottom Quark Production in pp 
Collisions at y/s = 1.8 TeV 
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Abstract 

This Letter describes a measurement of the muon cross section originat- 
ing from b quark decay in the forward rapidity range 2.4 < \y^\ < 3.2 in pp 
collisions at ^/s = 1.8 TeV. The data used in this analysis were collected by 
the D0 experiment at the Fermilab Tevatron. We find that NLO QCD cal- 
culations underestimate b quark production by a factor of four in the forward 
rapidity region. 
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Measurements of b quark production at the Tevatron have provided valuable information 
in the study of perturbative QCD. Cross sections measured by both the D0 [2,0, and CDF 
|| collaborations in the central rapidity region (\y b \ < 1.5) are systematically higher (by a 
factor of two to three) than the nominal values predicted by next-to-leading order (NLO) 
QCD j4|. Furthermore, discrepancies between some of the predicted and measured shapes 
of bb correlation distributions J5] indicate that the difference between data and QCD cannot 
be explained by a simple normalization factor. 

CDF has recently measured the bb cross section in which one quark is produced in 
the forward pseudorapidity region (1.8 < \r] b \ < 2.6) and the other in the central range 
(\r] b \ < 1.5) ||. This measurement is a factor of 2.4 higher than the NLO QCD prediction. 

Our measurement of the forward cross section of muons originating from b quark decay 
extends these studies to the previously unexplored rapidity region (2.4 < |y M | < 3.2), and 
provides further insights into the discrepancy between b quark production measurements 
and theoretical predictions. 

Forward muons are measured by the D0 detector [7] using the Small Angle MUon Spec- 
trometer (SAMUS) |S|||. SAMUS consists of two identical systems, each with three drift 
tube stations and a 1.8 T magnetized iron toroid, on either side of the interaction region. 
The momentum resolution of this system varies from « 19% at 20 GeV/c to « 25% at 
100 GeV/c. Muons reaching the SAMUS chambers traverse approximately 20 interaction 
lengths of material, reducing the hadronic punch-through background to a negligible level. 
This region does, however, face a large combinatoric background due to the flux of beam jet 
related particles. There are on average 6 to 14 hits per plane in a given bunch crossing, and 
the drift tubes near the beam axis have an approximate 5% occupancy. 

The data used in this analysis come from special runs taken at low instantaneous luminos- 
ity during the 1994-95 collider run. The integrated luminosity for these runs is 104 ±6 nb . 
The trigger required the presence of an inelastic collision near the center of the detector and 
at least one track in the SAMUS detector with an apparent pj- > 3 GeV/c pointing back 
to the interaction region. Muon candidates were also required to have an associated energy 
deposition in the calorimeter. The hit multiplicity in each layer was also required to fall 
below a maximum cutoff to improve background rejection and lower the trigger rates to an 
acceptable level. 

Muons are selected offline in the rapidity range 2.4 < \y^\ < 3.2, with p^ < 150 GeV/c 
and pj. > 2 GeV/c. Single interaction events are selected by requiring only one reconstructed 
vertex in an event, leaving an effective integrated luminosity of C = 75 ± 7 nb - . Muon 
tracks are required to have at least 15, out of an average of 18, hits. To ensure a good 
momentum measurement, we require muons to traverse a magnetic field integral of at least 
1.2 T-m. Muons are also required to be associated with a track-like object in the calorimeter 
with energy deposition consistent with that of a minimum ionizing particle. With these cuts, 
the combinatoric background is determined using both data and Monte Carlo (MC) to be 
less that 1%. The number of surviving muons in this sample is = 5106. 

The muon trigger and track reconstruction efficiencies are obtained using data and MC 
single muons, with detector simulation using GEANT |H| , superimposed onto real minimum 
bias events. The trigger efficiencies for the hit multiplicity cut [(31±2)%] and the calorimeter 
confirmation [(95 ±1)%] are obtained from data, as are the offline cut efficiencies for energy 
deposition [(94 ± 3)%] and number of hits on a track [(96 ± 2)%]. The overall detection 
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FIG. 1. The inclusive muon cross section in the forward region as a function of pj, (per unit 
rapidity). The dashed line shows the expected contributions from ir/K decays. 



efficiency is 1% for pj. = 2 GeV/c and reaches a plateau of 10% for pj. > 9 GeV/c. The 
MC momentum scale and resolution are shown to be correct to within 2% by comparing the 
peak values and widths of the reconstructed J/tp signal from data || and MC. 
The muon cross section is calculated as follows: 

df T dy» CAy^Ap^ e ' U 

where / smr is a correction factor that accounts for momentum smearing, and e is the detection 
efficiency. As there are high correlations between kinematic variables and cuts, / smr and e 
are determined by 

N"f ami 1 ff(data)ff(MCgen) 

e e data tf(MCreco) ' 1 j 

where edata is the combined data-based efficiency of the previously described cuts not simu- 
lated in the MC, and the H's are matrices with elements corresponding to two-dimensional 
histograms in the (pj^y^ 1 ) plane. if(data) is the data distribution after all offline cuts; 
-ff(MCgen) is the generated Monte Carlo distribution, and if (MCreco) is the reconstructed 
MC distribution with full detector simulation and the same cuts as the data. The histograms 
are segmented with 25 bins in pj> from to 25 GeV/c, and 7 bins in rapidity from 2.0 to 
3.4. The MC events are weighted in an iterative procedure to match the corrected pj. and 
rapidity distributions of the data. This method is found to give consistent results (within 
3%) regardless of the shape of the initial distribution. The resulting reconstructed MC dis- 
tributions also agree quite well with those of the data for all kinematic variables of interest 
after the weighting procedure. 
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TABLE I. Forward muon cross sections (per unit rapidity). 
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3-4 


3.4 


282.5 


± 


7.5 


± 


45 


92.2 ± 


33.1 
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0.078 


17.6 


± 1.0 
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0.675 ± 


0.081 


3.94 


± 0.36 


± 


0.83 


8-9 


8.4 


3.17 


± 


0.34 


± 


0.63 






0.685 ± 


0.075 


2.17 


± 0.23 


± 


0.50 
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0.067 
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0.057 ± 


0.011 ± 


0.026 






0.749 ± 


0.062 
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The inclusive muon cross section in the forward rapidity region (which includes both 
muon charges) is shown in Fig. [l| and Table [| The systematic errors in this measurement 
vary as a function of p T from 15 to 45%. They are dominated by uncertainties associated 
with the momentum smearing correction [(6 — 41)%], the single interaction luminosity (10%), 
and the trigger efficiency (8%). 

The contributions to this cross section from cosmic rays, hadronic punch-through, and 
W/Z decay are negligible (determined using both data and MC). The pion and kaon decay 



contribution is obtained using ISAJET [FJJ, which we find to be in agreement with the 



charged particle cross section measured in the central region ||12|| . The excess above the 
tt/K contribution is attributed to b and c quark decay. The fraction of this excess due to b 
quark decay (/&) can be obtained using the transverse momentum spectrum of the muons 
relative to that of an associated jet (p^ 1 ), but, because of our jet reconstruction threshold 
of Et > 10 GeV, only (7.9 ± 0.8)% of the events in the forward region have a reconstructed 
associated jet. We must, therefore, rely on a NLO QCD MC to determine /&. 

In this Monte Carlo, b and c quarks are generated according to the px and rapidity 



distributions of NLO QCD calculations J| using MRSR2 parton distribution functions JT3] 



quark masses = 4.75 GeV/c 2 and m c = 1.6 GeV/c 2 , with renormalization and factoriza- 



tion scales fi = /Jo = y^ 2 + Pt- The four-momenta of the quarks are input to an ISAJET 
MC which simulates initial and final state radiation, as well as quark fragmentation and 
decay. The theoretical uncertainty is determined by varying the parameters rrib from 4.5 to 
5.0 GeV/c 2 , m c from 1.3 to 1.9 GeV/c 2 , and /i from fi /2 to 2fi . The Peterson fragmen- 
tation parameters |TJ]]( e & = 0.006, e c = 0.06) are also varied by 50%, as are the branching 



ratios within their errors [15]. This simulation predicts that 8.5% of the muons should have 



a reconstructed associated jet, which is consistent within errors with what is found in the 
data. 

We check the validity of this MC by comparing its prediction for to that determined 
from our entire 1994-95 data set. Thirty-one thousand forward muons with an associated 
jet are selected from low pt single muon and muon+jet triggers. The trigger requirements 
keep the physics content of this sample the same as that of the cross section sample. The 
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FIG. 2. Data p^ 1 distributions for two selected p 1 ^ ranges. The solid line shows the fit to the 
data, with broken lines showing contributions from b quark (dashed), c quark (dotted), and tt/K 
(dot-dashed) decay. is the b quark fraction after tt/K subtraction (errors are statistical only). 

full sample is unsuitable for a cross section determination, however, as there is a large 
uncertainty in its normalization due to the various trigger thresholds and pre-scales, and 
luminosities that the data was taken with. The b quark fraction is determined by fitting 
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FIG. 3. fb for muons with an associated jet as measured from data p^ 1 fits (triangles and circle) 
and as predicted by the NLO QCD MC (dot-dashed curve) The prediction of /& for muons without 
the jet requirement is shown by the solid curve with uncertainties indicated by dotted curves. 

the p™ 1 distributions (in various ranges of pj) to the expected shapes from b quark, c quark, 
and tt/K decay (see Fig. as determined from ISAJET MC. The shape for tt/K decays 
was found to agree with the data distribution sample in the p^ range 0.5-1.0 GeV/c which 
is dominated by these decays. As is shown in Fig. |3], the NLO QCD Monte Carlo agrees 
quite well with the measured fb obtained in the p^ 1 fits of both the entire data sample, and 
the subset of events from the cross section sample that have a jet associated with a muon. 
Having shown that the MC is reliable for events with muons with jets, we assume it is also 
reliable for inclusive muons. 

Subtracting the tt/K contribution from the inclusive muon cross section and multiplying 
the result by the QCD MC predictions for fi, gives the cross section for muons originating 
from b quark decay. Our measurement, which includes both muon charges, and sequential 
b — > c — > n decays, is shown in Fig. |4] and Table |. 

The systematic uncertainties of this measurement include those of the inclusive muon 
cross section, with additional uncertainties due to /& and the tt/K subtraction. The con- 
tribution to the muon cross section from tt/K decay is predominantly in the low p£ bins. 
Conservatively assuming that the data in the 2-3 GeV/c bin (see Fig. [1]) is entirely due 
to tt/K decay, we determine that the ISAJET normalization is correct to within a factor of 
1.35. This factor is used to determine the uncertainty in the higher pj, bins. 

Also shown in the figure is a cross check of our measurement. We determine the cross 
section using the same events, but now require the muon to be associated with a jet, and use 



8 



10 



2 



> 

CD 

O 10 



2.4 < ly M l < 3.2 

a inclusive muon 
o muon + jet 



£Z 



a. 

Q. 



=i 10 



-1 



b 



10 



-2 




NLO QCD, MRSR2 
m b =4.75 GeV/c 2 

/ 2 2 \ 1 /2 



4 5 6 7 8 9 10 



20 



p T " (GeV/c) 



FIG. 4. The cross section for muons from b quark decay as a function of pj> (per unit rapidity) 
as measured with the inclusive muon sample (triangles) and its sub-sample of events that have 
a jet associated with the muon (circles). The solid curve is the NLO QCD prediction, with the 
dashed curves representing the theoretical uncertainties. 

the values for /& that were determined in the p™ 1 fits to the entire data sample. We obtain 
the same cross section (within statistical errors) as we do in the inclusive muon analysis. 

The NLO QCD predictions for the forward muon cross section from b quark decay are 
also shown in Fig. || as a function of They match the shape of the measured cross section 
fairly well, but are approximately a factor of four lower than the data. 

By combining the forward cross section with that of a previous D0 measurement in 
the central rapidity range (\y^\ < 0.8) fl|] we can study the rapidity dependence of b quark 
production. Our measurement of the cross section for muons from b quark decay as a function 
of rapidity (da^ /d\y^\) is shown in Fig. |^ for both pj. > 5 GeV/c and > 8 GeV/c. The 
ratios between data and theory are shown in Table ||. We find that next-to-leading order 
QCD calculations do not reproduce the measurements. 
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FIG. 5. The cross section of muons from b quark decay as a function of \y^\ for pj,>5 GeV/c, 
and pt^ > 8 GeV/c. The solid curves are the NLO QCD predictions, with uncertainty bands shown 
by the dashed lines. 



There have been some recent theoretical attempts to account for the discrepancy between 
data and theory. New calculations based on a variable flavor number scheme |16| predict an 
increase in the b quark cross section by a factor of 1.2 - 1.5 with respect to the standard 
calculations which use a fixed flavor number scheme. An increase in the 5-meson cross 
section of 50% in the forward region and 30% in the central region can also be obtained by 



using a stiffer b quark fragmentation function than the standard Peterson form [17 . Neither 



of these effects, however, can bring the predicted cross sections up to the measured values. 

In summary, we have measured the inclusive muon cross section, and the cross section 
for muons originating from b quark decay, in the forward rapidity region of 2.4 < \y^\ < 3.2. 
We find that next-to-leading order QCD calculations underestimate b quark production by 
a factor of four in this region. 
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TABLE II. The cross section of muons from b quark decay compared to NLO QCD. Errors are 
statistical and systematic added in quadrature. 



p£ > 5 GeV/c 


rapidity (y) 


measured 
< (nb) 


theory 
< (nb) 


ratio 


0.00 - 0.80 0.40 


89 ± 16 


36 


2.5 ± 0.4 


2.40 - 2.65 2.53 


43.5 ± 9.4 


12 


3.6 ± 0.8 


2.65 - 3.20 2.85 


30.5 ± 6.6 


8.4 


3.6 ± 0.8 


p£ > 8 GeV/c 


rapidity (y) 


measured 
< (nb) 


theory 
< (nb) 


ratio 


0.00 - 0.80 0.40 


20.1 ± 3.7 


6.6 


3.0 ± 0.6 


2.40 - 2.65 2.53 


7.9 ± 2.2 


1.6 


4.8 ± 1.3 


2.65 - 3.20 2.84 


4.1 ± 1.1 


0.99 


4.0 ± 1.1 
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